Universal Attractors of Reversible Aggregate-Reorganization Processes 
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We analyze a general class of reversible aggregate-reorganization processes. These processes are 
shown to exhibit globally attracting equilibrium distributions, which are universal, i.e. identical for 
large classes of models. Furthermore, the analysis implies that for studies of equilibrium properties 
of any such process, computationally expensive reorganization dynamics such as random walks can 
be replaced by more efficient, yet simpler methods. As a particular application, our results explain 
the recent observation of the formation of similar fractal aggregates from different initial structures 
by diffusive reorganization [Filoche and Sapoval, Phys. Rev. Lett. 85, 5118 (2000)]. 
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Many complex systems like microbial colonies, col- 
loidal aggregates, or particle- adsorbate structures are 
composed of a large number of similar units . For the 
growths of such particle clusters, the model of diffusion- 
limited aggregation (DLA) has been introduced ||2| and 
intensively investigated in several generaHzations [jj |, ||. 
In these non-equilibrium models, particles are perma- 
nently injected into the system and stick irreversibly to 
an existing cluster which grows towards a fractal struc- 
ture. In real systems, however, particles are often re- 
arranged, dynamically transforming the aggregate struc- 
ture. This rearrangement is appropriately modelled by 
processes exhibiting reversible dynamics 10, H, 

Recently, Filoche and Sapoval have proposed a specific 
model of diffusion-reorganized aggregation (DRA) jH) in 
which particles detach from the boundary of a given clus- 
ter of fixed size and reaggregate after a random walk on 
the underlying square lattice. In numerical simulations, 
they observe the formation of fractal structures similar 
to those found in various irreversible models. The DRA 
process transforms different initial aggregates, e.g. sim- 
ple one-dimensional and two-dimensional structures, into 
statistically similar fractal aggregates after many reorga- 
nization steps. 

In this Letter, we analyze a general class of reversible 
processes that dynamically reorganize a fixed number of 
particles constituting one connected aggregate on some 
lattice. Particles are allowed to reorganize by any re- 
versible process that retains a connected aggregate, in- 
cluding diffusion as a particular case. We prove that all 
these processes exhibit equilibrium distributions that are 
globally attracting. Our general mathematical perspec- 
tive reveals that such attractors are universal, i.e. inde- 
pendent of the dynamical model details: In discrete-time 
models, this distribution only depends on the mechanism 
of particle disaggregation but is independent of the sub- 
sequent particle repositioning dynamics. In particular, 
this implies that ordinary diffusion and different mech- 
anisms like surface diffusion ^ or Levy-flights |^ are 
equivalent in this context. Furthermore, we show that 
in continuous-time models equilibria are independent of 



both the mechanism of repositioning and that of disaggre- 
gation. On a d-dimensional cubic lattice, if all connected 
aggregates are accessible by the repositioning dynamics, 
this leads to the uniform distribution that is known as 
the ensemble of lattice animals (cf. |Q, ^, ^), currently 
also investigated in number theory under the name poly- 
ominoes (||l^ and refs. therein). 

In addition to these analytical results, we present 
strong numerical evidence that, for sufficiently large 
aggregates, there is no significant difference between 
equilibrium distributions obtained from discrete- or 
continuous-time models. Our results are vaHd for a gen- 
eral class of processes, including reorganization of aggre- 
gates consisting of different kinds of particles as well as 
repositioning of whole clusters instead of single particles. 
Here, for simplicity of presentation, we focus on processes 
in which single, identical particles are reorganized. 

We consider reorganization processes of aggregates - 
certain connected collections of a fixed number N of par- 
ticles on some lattice. Starting from an aggregate A, 
one step of the discrete-time dynamics is defined by the 
following two rules (Fig. |l|): 

(i) Randomly select a particle from a set of free par- 
ticles F{A) with equal probability, where F{A) is 
chosen such that the aggregate remains connected. 

(ii) Reposition this particle by a reversible process at 
some site x of possible aggregation sites D{A, B) of 
an intermediate object. By reversibility, this par- 
ticle is again a free particle of the new aggregate 
B. 

The above two rules are repeatedly appHed to the evolv- 
ing aggregate. Note that physically the step-by-step ap- 
plication of the two rules refiects a separation of time 
scales that is realized in many processes ||l||: The disag- 
gregation of a particle occurs much slower than its repo- 
sitioning. To model a specific reorganization process, the 
sets F{A) of free particles (i), the repositioning (ii) as 
well as the underlying lattice can be defined as desired 
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Figure 1: Model of a specific reversible aggregate- 
reorganization process: One particle (checkered) is selected 
from a set of free particles F(A) (gray) of aggregate A ac- 
cording to rule (i) and detaches from the aggregate. By a 
reversible process, e.g. diff^usion, it is repositioned at one (x) 
of the free aggregation sites D{A,B) (•) of the intermediate 
object according to rule (ii). Reversibility requires that the 
deposited particle (checkered) is again a free particle in F{B) 
(gray) and may be selected in the next reorganization step. 



as long as every aggregate under consideration can be 
reached iterating the above two rules. 

To be specific we focus on d-dimensional cubic lattices, 
but the analysis can be analogously applied for arbitrary 
lattices. We fix aggregates on a lattice of {N + 1)'^ sites 
with periodic boundary conditions and view translated 
aggregates of the same shape as different states. A given 
process is identified as a Markov chain 0| with state 
space Sn consisting of all connected aggregates of N par- 
ticles, which have at least one free particle, i.e. non-empty 
F{A). The state of this discrete-time process is then de- 
scribed by the distribution pt{A), giving the probabil- 
ity that after t reorganization steps the aggregate A is 
present. Starting with initial distribution po{A), a tran- 
sition matrix T specifies the time evolution 



Pt = Tpt-i 



(1) 



of the state vector pt := (pt(Ai), . . . ,pt{A\s„\))'^ ■ As Sn 
is finite, the long-time {t oo) behaviour of the discrete- 
time Markov chain is determined by a steady state dis- 
tribution p2 — Tp^. If the chain is irreducible, i.e. every 
aggregate A g Sn can be reached by the process from 
every other, p^J is the global attractor towards which ev- 
ery initial distribution converges This uniqueness 
property depends on the choice of free particles F{A). In 
general, one reasonably restricts the state space ^jv to an 
appropriate subset of aggregates consistent with F{A). 
For instance, if F{A) is the set of all one-bond parti- 
cles [||, only tree-like aggregates are considered. In such 
cases our analysis applies to the respective irreducible 
subset of Sn (also simply denoted by Sn in the follow- 
ing). With this convention, the stationary distribution 
p)5 is the unique attractor for every initial distribution 
Po on Sn ■ We now specify the transition matrix T which 
completely determines the dynamics of the process. 

If a one-step transition from aggregate A to B is pos- 
sible, denoted by ^ i?, the particle which is disaggre- 
gated and its reaggregation site are determined uniquely 



||l2| . Thus the underlying lattice does not explicitly en- 
ter the analysis. The transition probabilities are given 
by the matrix elements 



'BA 



\F(A}\ "^BA 





\{A ^ B 
otherwise 



(2) 



which, if non-zero, consist of two factors: According to 
rule (i) a free particle is randomly selected from F{A) 
with uniform probability 1/\F{A)\. Then, the particle is 
repositioned at the intermediate aggregate with probabil- 
ity tba to obtain aggregate B according to rule (ii). The 
reversibility of the repositioning procedure is character- 
ized by the symmetry tba = "tab- Thus, the dynamics 
is described by the master equation 



Pt+i{A) ~ pt{A) = [TABPt{B)-TBAPt{A)] 

Pt{B) pt{A) 



E 



\Fm \FiA)\ 



TAB .(3) 



A steady state solution of this equation is obviously given 
by the equilibrium distribution 



(4) 



BeSn 



which fulfills the condition of detailed balance, i.e. 



\F{B)\ 



\F{A)\ 



for all A,B £ Sn- Prom the above we 



know that p^ is the unique attractor of the dynamics 
(^. Hence, in equilibrium the probability to observe an 
aggregate A is proportional to its number of free 

particles. 

It seems natural to reaHze repositioning by a random 
walk of the disaggregated particle, because it models dif- 
fusion in real systems |^, |^. This random walk deter- 
mines the transition probabilities tab — tba between 
aggregates A and B. However, independent of the ac- 
tual values of tab , and hence the repositioning procedure 
(ii), the process leads to the same equilibrium distribu- 
tion p^. Thus, to study the equilibrium properties of the 
reorganization process, one can replace the random walk 
in numerical simulations: For instance, after selection of 
a particle, it may just be repositioned at every aggre- 
gation site with equal probability tab — ^/\F>{A, B)\, 
where D{A, B) is the set of possible aggregation sites 
of the intermediate aggregate (cf. Fig. || and ref. |0|). 
This straightforward procedure is much simpler to im- 
plement than random walks and significantly reduces the 
computational effort. Furthermore, one expects a faster 
convergence to the equilibrium distribution ||l^: While 
for diffusive repositioning the particle re-aggregates with 
high probability near or even at its disaggregation site, 
the alternative process leads to a fast spreading of parti- 
cles so that the state space is likely to be sampled faster. 
Moreover, diffusing particles in d > 2 dimensions have 
an infinite average return time leading to computational 
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complications ||5|, |g| that can be avoided using alternative 
processes. 

Up to now, we have considered discrete-time evolu- 
tion rules, which arise when performing computer sim- 
ulations. However, for real reorganization processes, 
continuous-time modelling, where the free particles F{A) 
disaggregate independently at a rate 7, often seems more 
appropriate. This leads to a continuous-time Markov 
chain with completely symmetric transition rates 'Jtab 
occurring in the master equation 

jPtiA)^ J2 7 [Pt{B)-pt{A)]rAB- (5) 

Thus, the attractor of these continuous-time processes is 

p:{A) = 1/\Sn\. (6) 

This is a universal equilibrium distribution which is in- 
dependent of all dynamical model details. It is important 
to note that, although F{A) does not enter (||), it im- 
plicitly determines the set Sn of accessible aggregates. 
In the case of d-dimensional cubic lattices and for every 
choice of F{A) for which the process can reach all pos- 
sible connected aggregates of N particles (cf. our simu- 
lations below), this uniform distribution p* is known as 
the ensemble of lattice animals |lO| . 

The question arises, how relevant the differences be- 
tween discrete- and continuous-time modelling are. To 
answer this question, we consider the fraction 

MA) = \F{A)\/N (7) 

of free particles of an aggregate A of N particles. 
This observable determines the difference in the equilib- 
rium distributions of discrete-time and their associated 
continuous-time models because it directly enters (Q) in 
the form p*aiA) = fN{A)/Y,BfN{B). In equilibrium, 
the average fractions of free particles obviously satisfy 
the inequality 

(/A.)p. - E fN[A)p*M) = > (/A^)p. (8) 

because (/Ar)p. — J2b fNiB)/\SN\ by definition. Never- 
theless, Monte-Carlo simulations of a specific process jl^ 
show that the averages (/at) clearly converge towards the 
same value for discrete- and continuous-time models (Fig. 

Differences become insignificant already for moder- 
ately sized aggregates. 

This is quantified by considering the widths of the equi- 

librium distributions a{fN) := ^(/w - (/jv))^^" (Fig- i, 

inset). We find a power law decay (j{fN) ^ N^" with 
a = 0.50 ± 0.02 for both time scales as would be ex- 
pected for large random structures. This has important 
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Figure 2: Comparison of average fractions (/jv) of free par- 
ticles: discrete-time (x) and its associated continuous-time 
model (O)- Error bars denote standard error of sample mean. 
Inset: Widths cr(/]v) for both time scales. 

consequences: Calculations similar to the above (^) show 
immediately that 

(5iv)p. = (3^)p. (1 + OiN-)) (9) 

for every observable gN{A). Thus, for sufficiently large 
aggregates, measurements of arbitrary observables, and, 
equivalently, the equilibrium distributions, in discrete- 
and continuous-time models cannot be distinguished. 
This refiects the fact that the number |i^(^)| of free par- 
ticles of almost all aggregates A G Sn is close to its 
uniform ensemble average {\F{A)\) for sufficiently large 
N. Since this fact is of purely configurational origin, we 
expect a similar convergence (||) for other models, also in 
d > 2 dimensions. 

As one important observable, we have studied the frac- 
tal dimensions of aggregates from both processes. Sample 
aggregates as well as box-counting dimensions (Dpf) are 
shown in Fig. ^. Even for moderately sized aggregates, 
no significant differences have been observed. With in- 
creasing aggregate sizes, the (Dn) approach from below 
Jl5[ the well-estimated fractal dimension D « 1.56 [Tc[| 
of square-lattice animals. 

The above results predict that the uniform distribution 
among all connected aggregates which can be reached 
by a specific process is the globally attracting, unique 
equilibrium distribution for a wide class of models. In 
particular, our results explain the recent observation 
that fractal aggregates are reached from arbitrary ini- 
tial structures by a certain diffusive reorganization pro- 
cess 1^. Nevertheless, our results strongly indicate that 
the equilibrium aggregates from that study are in the 
same universality class as lattice animals and thus have 
a fractal dimension D « 1.56 and not, as suggested, 
D = 1.74 ± 0.02 1^. A possible explanation of this dis- 
crepancy is that simulations in that study did not reach 
equilibrium at the time of measurement, despite a huge 
number of simulation steps. In fact, in our simulations 
Jl^ aggregates of = 10000 particles only reached equi- 
librium after about 2 x IG'^ particle rearrangements. 
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Figure 3: (a,b) Equilibrium fractal aggregates of A'^ = 10000 
particles, started from a 100 x 100 square, after 10^ particle 
rearrangements |^^ . (a) Discrete-time model, (b) associated 
continuous-time model, (c) Box-counting dimensions (Dm) 
for both time-scales are indistinguishable and approach from 
below the fractal dimension D ^ 1.56 (dashed line) ||9|, p^ . 



consisting of diflerent kinds a of particles disaggregating 
with particle-specific rates ja are also described because 
a transition A ^ B determines the disaggregating parti- 
cle uniquely. Moreover, the theory may be generaHzed to 
reorganization of whole sub-clusters rather than just sin- 
gle particles (cf. also |^ Hence more complex models 
of aggregate-reorganization can also be captured. 

We have imposed only a small number of requirements 
on the dynamics which appear to be natural assump- 
tions in many physical situations, in particular reversibil- 
ity of repositioning and separation of time-scales. Thus 
our main results might not only be of importance to 
aggregate-reorganization processes. We hope that the 
basic ideas can also find appHcations in studies of other 
complex systems exhibiting reversible reorganization. 

We thank M. Biehl, D. Brockmann, T. Geisel, L. Huf- 
nagel, T. Kottos, and M. Prahofer for useful discussions. 
This work has been supported by the Max Planck Soci- 
ety. 



It is important to note, that in certain systems such a 
slow convergence may imply that the equihbrium state is 
not reached during the time of real or numerical exper- 
iments (cf. 1^). Thus although an equihbrium distribu- 
tion is guaranteed to be reached, the long-term behavior 
of such systems may be governed by transient states with 
long lifetimes and non-universal statistical properties. 

In summary, we have shown that large classes of re- 
versible aggregate-reorganization processes possess uni- 
versal attractors which are independent of the model de- 
tails but may in general depend on whether discrete- or 
continuous-time modelling is used. Furthermore, our nu- 
merical investigations indicate that statistical properties 
are essentially identical for both time scales. These re- 
sults imply that commonly used diffusion modelling can 
be replaced by computationally fast methods which in 
addition lead to a faster convergence towards the (iden- 
tical) equilibrium distribution. In general, one specific 
process of aggregate-reorganization can be investigated 
in a number of different alternative models. Which alter- 
native is used can be freely chosen depending on which 
one is more convenient to study analytically or by com- 
puter simulations. 

In particular, our results explain the recently discov- 
ered "statistical attractors" in diffusive reorganization 
processes ||^ but predict a different final equilibrium 
state. More generally, the results can serve as a basis 
for investigating equilibrium statistical properties of a 
general class of aggregate-reorganization processes. Re- 
versible reorganization is always characterized by the 
symmetry tab = tba, which, in turn, has been the main 
requirement for universality. In particular, aggregates 
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